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ABSTRACT 


The operation of a bolometer thermal sensor is analogous to that of a charging | 
capacitor in a generic RC circuit. As such, circuits containing bolometers can be 
analyzed with standard circuit simulation programs such as PSPICE. This thesis deals 
with the development of a bolometer model by using PSPICE with the aid of Analogue 
Behavior Modeling (ABM) capability, which allows the user to program circuit 
components with basic mathematical functions. The predictions of the model were 
found to be in good agreement with the reported experimental data, which demonstrates 
the accuracy of the model. The model was used to design a self-heating compensated 
thermal sensor with enhanced signal integration capability to improve the signal-to-noise 
ratio. We believe the model can be used to analyze any circuit containing bolometers to 
optimize the performance. 
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I. INTRODUCTION 
Infrared (IR) technology remains vital to military, law enforcement and fire- 
fighting applications. Numerous means of detecting IR radiation exist, each with its own 
drawbacks and advantages. Among these, for example, the HgCdTe and InSb-based 
photodetectors typically have superior performance and moderate cost, but they, as well 
as many others, require cooling to reduce dark current (Rogalski, 1994). This in turn 


requires a larger, heavier overall package and higher power consumption. 


The need for small, lightweight, uncooled, sensitive IR devices has grown in 
particular in the military as a complement to the trend of smaller unit tactics and in 
response to the need for more intelligence and identification capability at the organic 
level. From imaging devices mounted on rifles to thermal detectors disguised as small 
plants providing perimeter surveillance, uncooled IR technology holds the promise of 


providing an edge over the enemy. 


The thermistor bolometer, which we will refer to simply as the bolometer, is a 
temperature-variable resistor which can be used in uncooled systems. By measuring the 
change in resistance, a signal can be found which is dependent only on the change in 
temperature which the bolometer experiences due to incident IR flux. A number of 
difficulties have been overcome, such as vacuum-encapsulating the bolometer to 
eliminate unwanted conduction of heat from the sensor by surrounding air. Nevertheless, 
a large obstacle to high sensitivity remains: the signal itself is affected by the very act of 


measuring it. 


To measure the change in resistance which has been caused by heating the 


bolometer, a bias voltage must be applied across it. The resulting current may then be 
l 








measured, and the change in resistance may be calculated. The difficulty in measuring 
the optically-induced change in resistance comes about because of the additional ohmic 


heating that the current generates. If we calculate the incident power due to a typical IR 


signal (~ nW) versus this self-heating (~mW), also called Joule heating, we find that at 


detector saturation the Joule heating is typically three orders of magnitude higher than the 


IR signal. 


The issue then becomes how to discriminate the signal from the self-heating 
response. In most applications, the self-heating response is first measured in the absence 
of an IR signal. This dark current value is then subtracted from subsequent measured 
signals. While this technique is relatively inexpensive and easily implemented, it is 


prone to error” and limits the dynamic range of the amplifying circuitry. 


A more promising technique is the use of a reference bolometer, which, if 
incorporated Aare a Wheatstone bridge arrangement (see figure II.7), can theoretically 
cancel out the self-heating response if it is optically isolated from incident IR radiation. 
One method was introduced in a European patent (Kimura, et al, 1997), which provides a 
physical shield for the reference bolometer. This technique may be effective, but it is 


very costly when implemented into micro-machined bolometer circuits. 


The necessity for low-cost, easily fabricated silicon chip technology prompted 
Karunasiri, et al, to investigate the possibility of rendering the reference bolometer 
translucent to IR by means of adjusting its thermal characteristics. Promising physical 


experiments validated this concept (Karunasiri, 1999), but it has not been fully explored. 


* an average value is subtracted from the actual integrated signal value, for example. 
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Analytical models of the process can be derived using first-order approximations, but 
what was needed next was a suitable computer model that could account for higher-order 
effects in order to further refine the concept of manipulating the reference bolometer’s 


thermal characteristics to improve self-heating cancellation. 


The purpose of this thesis is to develop a tool to conduct parametric studies of 
bolometer operation. What follows is the development of such a model using PSPICE, 
an industry-standard electrical engineering circuit analysis tool. PSPICE is used in a 
novel way to simulate the thermal characteristics of the bolometer. The model is verified 
against physical experimental data, and then an example is given of how the model can 
be used to explore and refine the concept of self-heating cancellation through the use of a 


compensating bolometer. 
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Il. THEORY OF BOLOMETER OPERATION 


A. GENERAL MODEL OF THERMAL DETECTORS 


1. The Heat Transfer Equation 


The operation of a generalized thermal detector can be understood by using a 
commonly accepted model wherein the sensor element is connected to a heat sink. 
Energy absorbed by the sensor is transported via a bridge with thermal conductance G to 
the heat sink. From the heat sink, the energy is either removed by a thermo-electrical 
(TE) cooling process or it dissipates as the heat sink re-establishes thermal equilibrium 
with its environment. 


sensor Heat Sink 








Figure IT. 1 Thermal detector model. 


A natural trade-off exists between the signal response level of the sensor and how 
quickly the sensor can acquire a new signal; both are dependent on characteristics of the 
sensor material as well as the physical connection to the heat sink. The process described 
here is governed by the heat transfer equation, which is given by (Boreman, Dereniak, 


1996) 


P Ser. (II.1) 


incident d 
t 
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The heat transfer equation, which is a statement of energy conservation, states that the 
incident power (Pinciden:) equals the thermal conductance (G) times the change in 
temperature, plus the heat capacity (#) times the rate of change of temperature. Thermal 
‘conductance, not surprisingly, is the ability of a device to conduct thermal energy and has 
units of watts per kelvin (W/K). Heat capacity (#1) is the ability of a device to accept 
thermal energy in the form of heat and has units of joules per kelvin (J/K). Typical 


values of G and H for a microbolometer are 10°’ W/K and 10° JK, respectively. 


Z Steady-state Operation 


Assuming a scenario wherein the sensor is subject to a constant incident radiative 
power flux, the sensor will reach a certain maximum temperature as its heat capacity is 
filled up. At this point the incident power simply passes through the sensor to the heat 
sink where it is dissipated. This change in temperature can be extracted from the heat 
transfer equation by setting the rate of temperature change equal to zero and solving for 


AT. 
When dAT/dt = 0, 


AT = (II.2) 


Q|'v 


3. Transient Operation 


Given the governing heat transfer differential equation H <* =P-—G-AT ,we 


wish to solve for the function AT in order to understand the bolometer’s change of 
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temperature as a function of time during the transient stages. Comparing the heat transfer 
equation to the equation governing the charging of a capacitor in a simple RC circuit with 
a constant current source J, note the similarity: 


yo 26 Kp=P versus CaS (11.3) 
dt dt R 


The RC circuit which equation (II.3) describes is shown below. 


Figure II.2 RC parallel circuit 
The solution to the RC circuit, assuming the initial value V=0 at t=0, is V = JR(1- e 7k ), 
and by comparison we should expect a similar solution to the heat transfer equation of the 


form AT = ail — e 74) if AT=0 at 0: 


Rearranging the heat transfer equation, 


dAT _ -dt 
GAT-P H- 


Let y= GAT -P and dy = GdAT . Then, 


%G nat 
y AL 


Moving G to the right side and integrating both sides, 











Sige 
Assigning Hy. = 7 , the time constant, 


| t : 
In(y) = ss +c, where c is an arbitrary constant. 


To solve for the constant of integration, we apply the boundary condition that at 
t = 0 the change in temperature (AZ) also equals zero, and hence y = -P (from 


y = GAT -P): 
In(-P)=c, andtherefore In(y)= a +In(—P). 
fy 


Collecting terms, 


and solving for y viele 
y= _Pet 
Substituting the original value of y = GAT-P we arrive at 
GAT -P=-Pe”*, 


and finally 


_Pl,_o% .H 
aT ==(I-¢ \, where T = Ty. (11.4) 


Typical values of t are on the order of 10° seconds” and are obviously material and 


design specific. 











If the incident radiation is removed after the sensor reaches saturation, P = 0, and 


we may observe how the device cools: 


0 =G-AT + H(dAT/dt). 


Rearranging, 

AAT _ ~dt 

GAT A 
and 

d 

% acces 

y HO 
substituting y=GAT this time, . 
P - 
AT =-ae", (II.5) 


where P is the power used to achieve the steady state from which the temperature 
then falls. The following figure illustrates the relationship between the change in 
temperature and the application of constant power which is then “tumed off” 


instantaneously (step response). 


exponential 
decay region 


exponential 


growth region steady state region 


Change in Temperature 





Elapsed Time 


Figure II.3 General thermal detector step response. 
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4, Practical Considerations and Trade-Offs 


As equation II.1 shows, a bolometer measures incident power by changing its 
temperature. The physical parameter that is actually measured is the temperature change 
AT. A value for P may then be extracted by applying equation (II.2), assuming G and H 
are known quantities. The other assumption made here is that we are measuring AT once 
the sensor has reached its steady state. Otherwise, as figure II.2 illustrates, we measure a 
time-dependent transient eipbige akin to equation II.4, which is a more complicated 


function of P. 


To describe how the sensor responds to incident IR, we define the temperature 
responsivity of the detector by 


ar 
a 
Noting that AT =P/G for steady-state conditions [equation (II.2)], we may re-write 


R, = (11.6) 


equation (II.6) as 


R,=— (11.7) 


and we see that responsivity is inversely proportional to G. 


Another important parameter for a detector is speed, or temporal frequency 


response. In the derivation of equation (II.4), we determined that the time constant 


= H ~ l e . . 
T= S/- , where t= Ms f and Af is the effective frequency response bandwidth. 


Considering this and equation (II.7) we see that speed and responsivity are inversely 
proportional. One gives up frequency response for greater responsivity in thermal 
detectors. Indeed, it is evident that the responsivity-bandwidth product R7Afis a constant 


fer a bolometer of a given heat capacity. 
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B. INTRODUCTION TO THE BOLOMETER 

1. Description 

A bolometer is a device whose electrical resistance changes as a function of its 
temperature. Bolometers are also known as thermistor bolometers, or thermally sensitive 
resistors. The bolometer’s ability to change resistance is described by its temperature 
coefficient of resistance (TCR), which will be hereafter denoted by a. Typical a values 
for metals material run in the region of approximately 0.3% change per degree celsius, 
but semiconductors may have values as high as 2% change per degree celsius. The value 
of o is used to determine resistance of a bolometer at a given temperature as: 


R(T)=R,(i+a-[T-h)). (11.8) 


where Ry is the bolometer’s resistance at ambient conditions, Tj is the ambient 
temperature and T is the instantaneous bolometer temperature. (7-Tp will be designated 
by AT hereafter.) 

A microbolometer is simply a miniature bolometer micro-machined or photo- 
etched onto a substrate, typically silicon. Figure [1.4 shows a pair of microbolometer 


images. 


1] 











Figure [1.4 Photographs of a low thermal conductance microbolometer (left) and a 
high thermal conductance microbolometer (right). (Ramakrishna, et al, 1999) 


The bolometers shown above are identical in every way with the notable exception being 
the addition of two wide connections to the substrate heat sink on the right bolometer. 
These additional bridges provide more material though which heat may flow to the heat 


sink. Hence, the right bolometer has a higher thermal conductance. 


2; Bolometer Design 


From both engineering and economical points of view, it is of great importance 
that the entire sensor be designed to fit onto a single silicon chip as an integrated circuit. 
Attempting to attach leads to individual circuit elements will otherwise dramatically 
reduce the manufacturing yield and increase the cost per detector. A cursory look at the 
micromachining process of a bolometer may provide a more thorough understanding of 


its physical characteristics. Figure [1.5 (below) depicts a typical design. 
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Anisotropic 
Etching 


Figure 11.5 | Microbolometer pixel design. 


The surface features are etched using standard lithographic processes. In order to 
thermally isolate the bolometer (with the exception of its bridges to the heat sink), the 
region below the bolometer is etched, for example, using tetramethyl ammonium 
hydroxide (TMAH) solution. This solution interacts with the p-type (100) silicon 
substrate in a useful fashion: it only etches the silicon in the 111 crystal plane direction. 
These 111 planes can be clearly seen in the following scanning electron microscope 
photograph (figure II.6). During measurement the device is placed under vacuum to 
eliminate any convective heat transfer. The importance of vacunim encapsulating the 
bolometer is demonstrated in figure II.7, which is a series of images taken by a focal 
plane array of bolometer pixels as it is placed under increasing vacuum. The image 
quality is improved because the bolometers suffer from increasingly less convection heat 


loss under increasing vacuum conditions. 
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Scanning electron microscope (SEM) image of a microbolometer showing 


Figure IT.6 


the etches along the 111 plane. (Gu, ez al, 1998) 
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C. READ-OUT CIRCUIT USING A WHEATSTONE BRIDGE 


1. Purpose 


The Wheatstone bridge (Dereniak, 1996) is commonly used for isolating circuit 
components and making measurements on them. As a first step in employing a 
bolometer to measure temperature or infrared radiation (IR), the Wheatstone bridge 
provides a convenient circuit architecture. As will be discussed in section C.7, using a 
second bolometer for compensation improves the output signal, but for now it is 


instructive to look at a Wheatstone bridge with a single bolometer element in it, as 


depicted below. 
Ve. 
Ro Ro 
V, V, 


R, R,(1 +aAT) 


Figure II.8 | Wheatstone bridge with one bolometer. In the text, Vpias will be denoted 
by V, and AVou will be denoted by Vou 


2. Signal Voltage as a Function of Temperature Change 


As previously stated, the bolometer’s resistance changes as a function of 
temperature in the following manner: R(T) = Ro(1+a-AT). Provided that no current flows 


through the nodes connecting V; and V2, we may consider both sides of the Wheatstone 
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bridge as voltage dividers with a bias voltage (V;) applied, in which case the voltages 


with respect to ground may be calculated as follows. 


V,=->y,, et, (11.9) 
2R, R+RO) 


In an actual circuit we may measure the voltage across the bridge (Vou), which can be 


calculated as 


. 


_Ro(l+aAT) — : (II.10) 


V =V.-V,=yV, 
oe Pees 2 


Vou =V, 


ul 


1+aQ@AT il 
2+QAT 2 








1+a@AT { 2-—aQAT ] 
Vu =V, ie 
2+QAAT\2-aAT) 2 


y up |2tQAT-@°AT* 1 
°|  4-@° AT? 2 


where the quadratic terms may be considered negligible if |@AT|«1. Then, | 


y. -v,] 2*a0r 
4. 2 


and finally we are left with a relatively simple expression for the signal voltage (V,,:) in 





terms of the change in temperature: 


V 1) 27). (11.11) 
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From equation (II.5) we understand that AT is an exponential in time with a linear 
dependence on incident power, which means that in the Wheatstone bridge, V.,; is also an 
exponential with a linear dependence on power. Two forms of power influence the 
temperature of the bolometer. Incident power is due to the incident IR radiation, which is 
what we wish to measure. The other form is dissipated power due to the bias voltage, 
which is applied in order to pass current through the bolometer in order to measure its 
changed resistance. 


3. Power dissipated across bolometer 


We now wish to consider more carefully how the heating of the bolometer occurs 
due to an external bias voltage (Vpias). This is known as Joule heating or self heating. To 
calculate the quiescent power dissipated across the bolometer, we square the voltage 
across the bolometer and divide by its resistance: 


V2 R(T 
Po. = Ren where V, ey, and R(T) =R,(1+@AT). 


RD) _y2 Ry(1+@AT) 72 


2°65 


5 OF RET I b = 2 
(R, + R(T)) R? +2R,(R,(i+@AT))+(R,(1+aAT)) 


Youle 4R2+4R2QAT + Ra? AT? ° 


Dividing through by Ro and treating the quadratic term as negligible for |@A7|« 4, 





~ _ Vy l+a@AT 


—- : 11.12 
Joule AR, 1+a@AT ( ) 
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and we see that the power dissipation due to Joule heating is independent of AT as long 


as |a@AT| <4: 


Vy? 


Proute = aR ‘ (1.13) 


Once steady state has been reached, all of the incident power will move through 


the bolometer to the heat sink, as can be described by the following set of equations: 





= = —GAT and ina steady state, oe =0,s0 


incident 
dt 


ie 


incident 


=GAT. (11.14) 


As previously stated, equation (11.13) expresses power dissipated by the 
bolometer due only to the bias voltage. In an application setting in which we wish to 
measure the incident IR radiation, we must recognize that the voltage response of the 
Wheatstone bridge will be due to the sum of the incident IR power and the dissipated 


power due to Joule heating. 


4. Wheatstone Bridge Response to Bias Voltage 


Before determining the effects of external radiation, it is instructive to look at the 
output voltage of the circuit as a function of time, due to an applied bias voltage alone. 


Figure IJ.8. Given equations (II.11) and (II.5), we may now write 


ae =“b9 Sit (16% | (11.15) 
4 G 
where ¢ is the time elapsed since the application of V,. In the absence of incident 


radiation, total power dissipated by the bolometer is given by the equation 
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Vy? 


ne IT.16 
Joule 4 R, ( ) 
and therefore the output voltage of the Wheatstone bridge is 
3 -f 
V 4 (1-¢ a (11.17) 
16-R,:G 


The response, not surprisingly, takes on the same graphical form as that of the 
exponential rise in temperature depicted in figure 113. Let us look more closely by 
inserting some typical values into equation (II.17). Figure IL8 represents the output 


voltage of just such a circuit (parameter values are given in the figure caption). 


bias voltage = 0.8 V 


Output Voltage (mV) 





0 100 200 300 400 500 600 700 800 900 1000 
Time (ms) 


Figure 11.9 | Output voltage for a one-bolometer Wheatstone bridge with 
G=3.2x107W/K, H=4x10°J/K, R=16kQ, and a=0.27%. 


For t>>1 (steady state) the output voltage asymptotically approached the value of the 


prefix term in equation (II.17), which is called the saturation voltage 


3 
Ve -a 


Venton oo II.18 
Saturation l 6 ; R, ; G ( ) 
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and in the example of figure II.9 1s calculated to be approximately 10.4 mV. Equation 


(II.18) may be used to determine the bolometer detector system’s G from a measurement 
of V saturation: 


It is also useful to determine an empirical value for the system’s H parameter. For 


t << T, we expand the exponential term to the first order and simplify: 


3 
Pee 1-(1-4) 
16-R,:G T 


3 
Y= Nya t 
16-R,-Grt 
and using equation II.4, 


3 
Ve 0% 


Vou =. (II.19) 
16-R,-H 


From this equation we observe a linear response for small ¢ which is independent of G. If 
we were to collect a number of data points for this time region, we would obtain a plot of 


the form found in figure II.10. 


Output Voitage (mV) 


Time (ms) 


Figure II.10 Depiction of the output voltage response for Wheatstone bridge over a 
period of time t«r. 
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The fitted slope of the corresponds to the coefficient of ¢ in equation (II.19), which can 
then be used to solve for H. Hence, one can determine both H and G with one set of 
measurements. 


5. Radiative Thermal Conductivity 


One important consideration should be noted, however, concerning 
thermal conductance. Heat transfer occurs in three forms: convection, conduction and 
radiation. With the bolometer under vacuum, convection is negligible. Nevertheless, 
while heat put into the bolometer will be conducted through the bridges attaching it to the 
heat sink, it may also radiate. When an object heats up, it radiates in accordance with the 
Stefan-Boltzman Law (Dereniak, 1996), which gives the radiant power as 


P 


rad 


= Aso(T’ -T,'), (11.20) 
where 7p is the ambient temperature, A is the radiant surface area, emissivity € 1s a 
function of the bolometer material, and the Stefan-Boltzman constant 
o =5.67x10° watt/(m?-K*). <A true blackbody has an emissivity of exactly one. 
Considering that most metals are more reflective than emissive, one would assume our 
titanium bolometer would have a relatively low emissivity. Because the bolometer is 
typically biased with low voltages, it seems unlikely that emissivity would play a role at 
correspondingly small temperature changes (e.g. less than 100 K), but it certainly 


warrants analysis. 


The question is whether or not enough power is radiated away from the bolometer 
(P,adiated) SO aS to make a significant impact on the measured signal. To help describe the 
radiative transfer of energy, we define thermal conductivity due to radiation Gyag as 
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P.. 
G ae radiated ; 9 l 
r ad A T ( ) 


with the understanding that it is not an actual physical parameter but rather a construct 
which can assist us in describing radiation in convenient terms. Substituting equation 
(11.20) for the radiant power into equation (2.21), 
G_W AT = Aco(T* -f,). (11.22) 
Substituting 79+AT for T, and moving a AT to the night side, 
Aeo([T,+AT] -1] 


rad A T 


Assuming AT is small relative to 7”, all but the first order terms can be dropped in a 


(11.23) 


4 id 


Taylor expansion of |Z, + AT] 








4 
[TZ + AT] = ty [ + a (11.24) 
0 
After the expansion and dropping of higher-order terms, 
[%+AT] = [ Fae | (11.25) 
Ty 
This approximation can be re-inserted into equation (II.23) to yield 
Aso G [ + — -1] 
6. | 11.26 
which results in 
G4. =4AeoT, . (11.27) 


This result surprised the author, but upon reflection it certainly makes sense. Every 
object at room temperature is radiating in accordance with the Stefan-Boltzman Law; a 
slight change in temperature will not have a major effect on that radiation. The fact that 


* At low bias voltages, typical AT’s are around 30 °K, compared with room temperature of 300 °K. 
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with small temperature changes G,aq is constant is fortunate for those who wish to 
establish a model that includes it. The method described earlier for calculating G remains 
valid, with the understanding that G = G,,,, sretive + Gradiative 

6. Figures of Merit 

The figures of merit for a general thermal detector discussed in chapter ILA.1 
may now be revised so as to describe the bolometer in a Wheatstone bridge circuit. We 


now define voltage responsivity as 


y 


R, =—*— (11.28) 
incident 
and from equation (II.2) we have 
AT se elon 
G 


for a saturated bolometer, which is what would be measured. Substituting equation (II.2) 


into equation (II.11), we have 


Vi, ak 





incident 11.29 
out 4 . G ( ) 
and equation (II.28) becomes 
Va Ty 
Ry = | py | | (11.30) 


Regarding speed, the same relationship holds true for the time constant t, which, as given 
by equation (II.4), is 7 = Hy. 

7. Measuring Incident Infrared Radiation 

In practice, if we were using our single bolometer Wheatstone bridge circuit to 
make measurements, we would expose it to the IR radiation for a more than three time 


constants in order to achieve steady state, apply a bias voltage pulse and measure the 


output voltage. If we compared this output voltage data with that of an experiment in 
23 














which we just measured the output voltage without incident IR, we would see something 


like the figure given below. 
Signal due to Self Heating and Incident oa 
& 
= Signal due to Self Haaeg 
> 
A 
6 





Time 


Figure JI.11 Exaggerated illustration of the output voltage to due Joule heating alone 
and to Joule heating combined with incident. IR. 


Figure IJ.11 illustrates the resultant relative output voltages and is only provided for 
pedagogical use. It is in fact a gross exaggeration; the difference between the two curves 
would not be discernable if the two curves could be drawn to scale. Typical values for 
incident IR fluxes are in the nanowatt (nW) range, compared to milliwatts (mW) due to 
Joule heating. In an actual circuit, the minuscule IR signal would have to be extracted 
from the Joule heating response prior to amplification. If it were not removed, 
introduction of the Joule heating signal into the amplifier would cause quick amplifier 


saturation. 


One way to compensate is to electronically subtract a calculated value of the Joule 
heating voltage from the composite signal. This method is used with numerous current 
devices, but it relies on analytically derived values to be subtracted from the experimental 
composite IR signal and Joule heating signal. If the signal is integrated for a certain 


period of time in order to improve the S/N ratio, cancellation becomes even more 
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difficult since the self-heating signal increases with time (see equation II.20) until the 
steady state. Another way would be to physically cancel the unwanted voltage by using a 
compensating bolometer on the opposite side of the Wheatstone bridge, as illustrated in 


the following schematic: 


circuit 





Figure II.12 Compensating bolometer Wheatstone bridge circuit. 


In this arrangement either bolometer B; or bolometer Bz is the “compensating” 
bolometer, which would mean that it should be optically masked from the incident IR 
radiation, but allowed to undergo normal Joule heating. With the objective of creating a 
micro-machined, integrated circuit on a silicon substrate, physically shielding the 
compensating bolometer from the IR radiation becomes impractical. However, by 
altering the thermal characteristics of the compensating bolometer, it can be made to 
appear virtually translucent to the IR irradiation. Specifically, by increasing its thermal 
conductivity, i.e. making the legs which saach the bolometer face to the heat sink larger, 
the energy which is absorbed will pass through the bolometer’s legs to the heat sink at a 


higher rate. 
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8. Compensating Bolometer Operations 


The use of a compensating bolometer with higher thermal conductivity, as 
opposed to an identical bolometer which is physically shielded from incident IR flux, 
introduces a host of dynamic changes to the operation of the Wheatstone bridge. 
Changes in signal strength as well as thermal time constants accompany any change in 
the ‘henna conductivity of the reference bolometer. Figure II.13 depicts the temperature 
changes experienced by a sensor bolometer B, (upper curve) and a reference bolometer 
B2 (lower curve) as they are both exposed to incident IR flux and then to a bias voltage 


step function. 


Change in Temperature 





Figure II.13 Temperature responses of a sensor bolometer and reference bolometer. 
The upper curve represents B; and the lower curve represents B2, where G2=2G}. 


In the interest of simplicity, incident IR power has been made equal to the power 


dissipated due to bias voltage. The thermal conductivity of the reference bolometer has 
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been arbitrarily selected to be twice that of the sensor bolometer, in order to demonstrate 
how thermal conductivity affects both the maximum signal strength and the speed with 


which the bolometer heats up, in accordance with equations (II.2) and (11.4), respectively: 


AT =—, (11.2) 


-P(,_.% _H 
aT ==(I-¢ where T = Un. (11.4) 


The signal we wish to measure is the maximum temperature of the sensor 
bolometer due to incident IR flux, i.e., the upper curve before bias voltage is applied. 
With the reference bolometer in place, the actual IR signal available is the difference 
between the two curves before the bias voltage (A7). From this perspective, the reference 
bolometer has degraded the signal (by half). In practice, if we use a reference bolometer 
with a thermal conductivity one-tenth that of the sensor bolometer, we will experience 
only a 10% signal loss. Hence, the first of two design constraints is evident: G2 should be 


as large as possible to make B2 optically transparent to IR radiation. 


The benefit of the reference bolometer becomes manifest when the bias pulse is 
applied and the two temperature levels rise exponentially. Considering the more realistic 
case in which ‘esis heating power is orders of magnitude higher than IR power, one 
might be tempted to wonder how an accurate signal can be extracted at all with the use of 
a reference bolometer. Analytic results in the case where G2=2G; and dissipated power 
equals 1000 times incident IR power are given in the following figure, generated by 


MATLAB. 


Zi 
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Figure II.14 Scale view of response in which applied power = 1000xIR power and 
G2=2xG;. The right hand figure illustrates the initial equivalent linear growth for both 
curves. 


The original signal is completely buried shortly after bias voltage is applied for the read-’ 
out. F ortunately, while the temperatures of B; and B2 grow with different exponential 
rates, there is a region in which their growth rates may be approximated as linear and 
independent of thermal conductance, which is illustrated in the right hand graph of Figure 
II.14. Rearranging the heat transfer equation, | 


d(ar)=—at-<ar. (11.31) 


If dt<<r, then equation (II.31) can be approximated as 


d(AT) = = dt (11.32) 


Therefore, in the region dt<t, where dt is the bias voltage pulse width, we should be 
able to integrate the signal to measure AT quite accurately, as long as the heat 


capacitances of the two bolometers are the same. Here we see our second design 
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constraint: for t to be large compared to dt, Gz must be small, so that it may heat up at 
the same rate as B; for at least long enough to obtain a read-out from integration. 

These two competing design constraints were the motivation for a computer 
model that could be used to theoretically optimize various Wheatstone bridge 
configurations. With this in mind, Section III deals with developing and verifying an 


acceptable model using an industry-standard electrical engineering tool called PSPICE. 
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Ill. MODELING THE BOLOMETER WITH PSPICE 


A. DISCUSSION OF PSPICE 


PSPICE is a computer simulation program that models the behavior of an 
electrical circuit containing analog devices. It is essentially a computer-based breadboard 
which allows one to predict AC and DC steady state and transient analyses. By using it 
to model the circuits described in Section II and then verifying the model against 
experimental data, we may confidently conduct parametric studies of the bolometer. 


1. Bolometer Heat Transfer Equation versus the Parallel RC Circuit 


As discussed in Section I, the heat transfer equation is of the same form as the 
differential equation for the parallel RC-circuit under constant current charging (Figure 


11.3.) The following are the two analogous equations: 


(III.1) 





Their analogous parameters are as follows, with their SI units in square brackets: 


Hk] >ClH} 
o[%|> 3 Vob 


rie] > V0) PL i 
P(t) [7] > 1(+) [4 


Using this analogy, we can design a circuit using PSPICE to mimic the temperature- 


dependency of the bolometer’s resistance. We can then insert this virtual circuit into 
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another circuit representing the actual Wheatstone bridge with a bias voltage applied to it. 
Special consideration must be given not to confuse the actual electric circuit values with 
the virtual “temperature circuit” values. The bolometer is modeled as an RC circuit with 
the resistor representing the thermal conductance and the capacitor representing heat 


capacity, as shown below. 





y2 
Fr ~ ome ) bolometer ~ uy Hy 


Figure II.1 ‘Virtual circuit representing a bolometer. 


Subscripts “R® and “V” stand for real and virtual quantities, respectively. R(T) bolometer iS 
given by R(T) = Ro(1+a-AT), and Voolometer is the value of the voltage drop across the 
bolometer in the Wheatstone bridge. The value of the current source in the virtual circuit 
corresponds to the power dissipated in the bolometer. In figure III.1, a voltage source 
maintains a virtual offset voltage of 300 V, placing the sensor’s initial temperature at 
room temperature conditions of 300 K. What remains is to connect the bolometer to a 


Wheatstone bridge and apply a bias voltage, all of which we can do using PSPICE. 
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B. THE SINGLE BOLOMETER WHEATSTONE BRIDGE 
i PSPICE Model of Single Bolometer Wheatstone Bridge 


The point of interface for the circuit described above is the bolometer itself, 
represented by an analogue behavior modeling (ABM) box, which allows the user to 
control the output through the use of the algorithms he chooses. As there is no 
temperature-dependent resistor in the ABM options provided by PSPICE, we must 
replace the function of the resistor with an available option: a variable current source, 
which we can program to set the current through to be the same as that through the 


bolometer it replaces. 






V(1)*V(1)/(Ro*(1400*(V(2)-V(3)))) 


= V(1)/(Ro*(1+0*(V(2)-V(3)))) = 
Figure III.2 | PSPICE model of Wheatstone bridge with a single bolometer. 


Figure III.2 represents the circuit as seen in the PSPICE format, with the exception being 
that the ABM formula accept only numeric values, as opposed to the formulas shown in 


the figure, which include a and Rp to help clarify the figure. The ABM box receives 


temperature inputs ([V (2)-V (3) | —> AT ) from the virtual circuit to insert the correct 
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current into the Wheatstone bridge, acting as if there were a temperature-dependent 


resistor there. The current generated is given by | 
V, 
Ttt)= ° Solometer ; 
) “Zp (7) (III.3) 


where R(T) is 


R(T)=R,(1+aAT). (14) 


The virtual circuit is able to calculate the temperature difference AT by inserting the 
voltage drop across the bolometer in the actual circuit into an ABM box (far left) to 
generate a virtual power, whose RC-circuit analogue, according to equation (III.2), is 


current. 


At this point the circuit has been designed to represent just the Joule heating effect 
without incident radiation. To include irradiance, a current source (analogous to power) 
would be placed in parallel with the resistor and capacitor in the bolometer circuit. 


2. Verification of the PSPICE Model and Early Comparison with the 
Analytic Model 


After creating the PSPICE model shown in figure III.2, the first step to 
verification is to observe the correlation of simulated data to the analytical predictions of 
Section IT. One of the first experiments run was to verify equation (II.13), which states 
that the dissipated power in the bolometer is constant, given a steady bias voltage (V;): 


2 
P V, 


=—., IJ.13 
Joule AR, ( ) 
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Recalling that a quadratic term of temperature (Rja°AT’) was dropped from the 


denominator, what this simulation really seeks to determine is for what region equation 


(IL.13) is valid. Using a bolometer with the same values as before (G =3.2x10"W/K, 
H=4x10°J/K, R=16kQ, and o=0.0027), the results for applied bias voltages 


ranging from 1.0 to 10.0 volts are given below. 





1.0 
Analytic values 
BEN cassette lest Sepals DataSet ts cette 
s 
x@) 
2 10.0 V 
© 
& 9.0V 
wm 04 8.0V 
ra 7.0V 
® 
3 6.0 V 
Oo. 0.2 5.0V 
4.0V 
3.0V 
2.0 V 
0.0 1.0V 
0 200 400 600 800 1000 


Time (ms) 


Figure III.3 Joule heating power dissipated for varying bias voltages. 


Also shown are the corresponding analytical values for the top three bias voltages (8.0, 
9.0, 10.0 V) as given by equation (II.13). What is quite apparent is the departure from the 
expected linear relationship. The nonlinear effects of the disregarded quadratic and 
higher temperature terms are in fact important. So what is the acceptable region of 


linearity? 
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We can determine the condition for acceptably disregarding the quadratic term in 
the equation 


R,(1+@AT) , 


oe eM anal ee 111.5 
4R?+4R-aAT+R2a7AT’ ° ae! 


by setting up the condition to be 


4RiaAT > Ra’ AT’ 


which then leads to 


is > AT. (II1.6) 
a 


In the case of a=0.0027 for Ti bolometers, the acceptable region must be AT «1480K. 


The exact region obviously depends upon how accurate a prediction one would like to 


make. The deterioration in linearity can clearly be seen in the following set of graphs. 


Power/Temperature Relationship for bias Voltage = 1.0V Power/Temperature Relationship for bias Voltage = 2.0V 
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Power/Temperature Relationship for bias Voltage = 5.0 V Power/Temperature Relationship for bias Voltage = 6.0 V 
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Power/Temperature Relationship for bias Voltage = 7.0 V ; Power/Temperature Relationship for bias Voltage = 8.0 V 
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Power/Temperature Relationship for bias Voltage = 9.0 V Power/Temperature Relationship for bias Voltage = 10.0 V 
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Figure IlI.4 Effects of increasing temperature on dissipated power. 
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One can actually begin to see degradation on the second graph (V;=2.0 V), and by the 
fourth graph (V,=4.0 V) there is significant deviation from the expected linear behavior. 
At V;=10.0 V, the steady state power is slightly over half of the value predicted by the 
analytic model. By this point the temperature change is at approximately 1600 K. Hence 
we have clear verification that, in order for our model to be accurate, the bolometer must 
be operated in the region well below 1480 K, as equation (III.6) states. 

5, Verification of the Model by Comparison with Physical Data 

Having established the acceptable temperature regime, we can also verify the 
model by comparison with actual physical data (Ramakrishna, et al, 1999). Two 
bolometers with distinct thermal properties were used in a ecm experiment to 
demonstrate the concept of a compensating bolometer Wheatstone bridge. The thermal 


characteristics are given in the following table. 





Bolomter B; Bolomter B> 
Resistance (at 294 K) 26kQ 94kQ 
a (at 294 K) 0.0027/K 0.0029/K 


Thermal Conductance (G) 3.61x10’ W/K 2.37x10° WK. 


Thermal Mass (H) 4.0x10° JK  —4.74x10* JK 
Thermal Time Constant (t) 111 ms 15 ms 
Responsivity (at 3 V) 6000 V/W 770 V/W 


Table III.1 Bolometer characteristics 


Below is a graph showing measured output voltage versus time for the two 
bolometers. Note that the upper and lower x-axes separately represent the two 


bolometers. 
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Figure III.5 | Experimental output voltages for the two bolometers as a function of time 
at 294 K (Ramakrishna, ef a/). 


Running a PSICE simulation with the same parameters, the following data is generated: 
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Figure III.6 Simulated output voltage for bolometer B, and bolometer Bz at 294 °K. 


The agreement between Ramakrishna and Karunasiri’s physical experiment and the 
simulated experiment of figure III.6 falls within 1%. The data for figure III.6 is collected 
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in PSPICE by placing “voltage differential probes” on either side of the Wheatstone 
bridge, in each case between the two resistors, as shown in figure III.2. We can simulate 
monitoring the temperature by placing the voltage differential probes on the upper and 
lower branches of the bolometer’s virtual temperature circuit in figure III.2. Doing so 


yields the following data. 


Time axis for B, (ms) 
0 20 40 60 80 * 100 


= V,= 0.8 V (Bolometer, ) 











poor V,= 1.3V (Bolometer,) 
Oo 


e) 


Temperature Change (K) 


0 200 400 600 800 1000 
Time axis for B, (ms) 


Figure III.7 | PSPICE simulated temperature change for bolometer B, 


While it is not possible to measure the steady-state temperature change in the physical 
experiment, it can at least be calculated analytically to determine how closely the model 
agrees with analytic predictions. For bolometer B;, with V,,;= 0.0091 V, V,=0.8 V and 
a = 0.0027 °K" , equation (11.11) yields AT = 16.85 K. The PSPICE simulation yielded a 
maximum temperature change of 17.02 K, which indicates agreement within 1%. For 


bolometer Bo, with Vou= 0.007 V, V>=1.3 V and @ = 0.0029, equation (II.11) yields AT= 
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Output Voltage (mV) 


7.427 K. In this case the PSPICE simulation generated a steady-state temperature change 
of 7.416 K, well under a 1% discrepancy. It should be noted that these slight errors may 
be due to estimating the saturation voltages from figure III.6, as actual numerical data 


was not available from the physical experiment. 


Before completing the compensating bolometer Wheatstone bridge, we can 
compare modeled data with the physical data on how the individual bolometers respond 


to Joule heating and incident IR in the dt <7 region, where dt is the bias pulse width. 
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Figure III.8 Physical data for bolometer Bj (left) and bolometer By (right). 


Once again, using bolometers with identical thermal characteristics, the following output 
voltage data is generated by PSPICE. (next page) 
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Figure III.9 PSPICE simulation of bolometer response to incident IR in addition to 
Joule heating. Bolometer B, is represented by the left graph, bolometer B by the right 


graph. 


In both figures JII.8 and IJI.9, we can see the IR signal riding on top of the Joule 
heating signal. Additionally, the IR signal is significantly stronger on bolometer, than on 
bolometer2, which was the goal of attempting to make bolometer, translucent to IR by 
increasing its thermal conductance. Nevertheless, because we are operating in the 
dt<<t region, the output voltages due to Joule heating for both bolometers are roughly 
equivalent, following the theory leading to equation (11.32). The PSPICE model appears 


quite accurate even at small transient times. 


4. Modeling the Compensating Bolometer Wheatstone Bridge and Final 
Verification. 


With the pieces in place and functioning properly, it remains to assemble the 


Wheatstone bridge with a compensating bolometer in place in order to cancel out the 
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majority of the Joule heating signal, leaving the IR signal and a roughly equivalent Joule 
heating signal as the input to the read-out circuitry, which would ostensibly begin with an 
amplifier (the very reason why it is advantageous to cancel out the enormous Joule 
heating signal). The following physical data show sood cancellation of the Joule heating 


signal. 


Without Compensation 


Output Voltage (mV) 


o = 280 nW 


@ = 110 nW 





0.0 0.5 1.0 1.5 2.0 


Time (ms) 


Figure IIJ.10 IR and Joule heating response of the bridge consisting of two bolometers 
contrasted with response of the bolometer B; at bias voltage of 3V. The two bolometers 
were exposed to the same amount of infrared from a blackbody source. 


wT 


The actual numerical physical data was not available, nor were the specifics of the 
experiment setup. For the PSPICE simulation, the Wheatstone bridge was set up with B, 


on one side with an equivalent resistor and Bz on the other side with an equivalent 
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resistor. (This setup is shown in figure JJI.12.) The simulation yielded the following 


data. 


o® 
2 ss 
10 RAO 
ad y” AY 
Without Compensation _.* ay? 
ov gee 
8 yo 2" Pra 
v a” 4 
Ss Oi " ao 
S tr aut 
E gy” 2. a° 
® 6 v a" ana 
mo y” . a> 
93) v” e* 2° 
5 oe? s* ao 
> 4 ore ta 
3 Os ." 4 - 
£ v” at ae , 
5 or” et ge? With Compensation 
2 ae so” 6° saa Side intents o = 280 nw 
Bc PP HI HHI AAAI EOS b = 440 AW 
: eoeececsessesooeee 


passbsdbdsdsbbbasd b =0 nw 


fErstorenceensnnsonssnveeereores 
4a 


YY YVYVYY VV YY TTL asssssbbdda 





Time (ms) 


Figure IJJ.11 PSPICE simulation of the IR and Joule heating response of the bridge 
consisting of two bolometers contrasted with response of the bolometer B, at bias voltage 
of 3V. The two bolometers were exposed to the same amount of simulated infrared flux 
from a blackbody source. 


The correlation between the physical data and the simulation data is favorable, 
once again under 1% difference. Nevertheless, while the comparison of the physical 
cancellation with the simulated cancellation are inconclusive, the individual bolometers 
matched very well, and it is assumed that combining the two onto one compensated 
Wheatstone bridge should not alter their individual performance, as each is subject to the 
same voltage. Therefore, on the basis of the excellent correlation between reality and 


simulation as demonstrated by figures ITI.8 and III.9, the model is deemed a success. 
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Figure III.12 Compensated bolometer Wheatstone bridge 


[Note: to include variable thermal radiation for cases when the condition AT «T, does not apply, ABM 
boxes should be generated in parallel in both bolometer circuits and programmed with the Stefan-Boltzman 
Law, equation (II.20).] 
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C. PARAMETRIC STUDY USING THE WHEATSTONE BRIDGE MODEL 


The virtue of the PSPICE model can now be realized in that it is readily possible 
to conduct research and design studies of microbolometers without actually 
manufacturing them, saving time and money. Given a certain initiadiin the designer 
may now quickly produce the necessary thermal characteristics required to design a 
useable device. It should be cautioned, however, that the assumptions and | 
approximations made in designing the model must still hold true in order to create an 
accurate design. What follows is an example of how bolometer applications can be 
improved by employing the compensation technique, as illustrated through the use of the 
PSPICE model. The parameter to be used for comparison will be integration time, which 
therefore warrants a brief discussion. 


1. Integration Time versus Noise Bandwidth 


A cursory look at the forms of noise affecting bolometers reveals that they are 
dependent on frequency bandwidth Af 
V,, =4kTR-Af Johnson noise 
Ve =i R? (ail) Vy noise 
where R is resistance, k is Boltzman’s constant, T is temperature, and B is a constant for 
My noise. (Shot noise is deemed negligible here). In our application, Afis determined 
by the pulse width of the sampling pulse, where pulse width directly correlates to 
integration time fint. Recalling from Fourier analysis, a narrower square pulse requires the 


addition of higher frequencies for its construction than a broader square pulse. Generally, 


the bandwidth is given by 
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2+ bint ) 





Af = (111.7) 


Hence, the longer one can integrate a signal, the less noise will remain. In the 
comparison which follows, the decision is arbitrarily made that the signal will be 
integrated for a time period sufficient for the voltage response due to Joule heating to 
equal the voltage response due to an arbitrarily chosen IR radiation flux. This concept 
reflects the feasibility that the subsequent electronic subeaeon of a Joule heating voltage 
on the same order of magnitude as the signal voltage reduces the probability of distorting 
the remaining signal voltage. 


Ze Application: Thermal Imaging Device 


The necessity and usefulness of lightweight, portable thermal imaging devices are 
obvious to military, law enforcement and fire fighting personnel. The ability to image 
targets, suspects, and fires or hot spots in completely dark and/or smoky scenarios is 
critical to the success of many kinds of operations carried out by such personnel. As 
explained in the introduction, there are currently numerous devices available, but 
typically the operator must accept a trade-off in terms of sensitivity versus weight, 
because the more sensitive devices typically need to be be cooled, which requires more 
power, weight and volume. The microbolometer fills this niche of high sensitivity and 
low weight/space at the expense of speed. Bolometers are relatively slow devices 
compared sith many others, but this lack of speed really only precludes them from high- 


speed applications, such as IR warhead seekers. 


Raytheon Systems Company has developed a prototype infrared imaging rifle- 
sight using an uncooled microbolometer focal plane array (FPA) (Radford, et al, 1998). 


(Implemenatation of the bolometer into a FPA is introduced in Appendix A.) This 
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system is typical for the kind of application needed in the Department of Defense. 
Because it is an integrated, micromachined chip technology, it is relatively low in cost 
compared with other microbolometer applications that require complicated 
manufacturing procedures to hide the reference bolometer from IR radiation, for 
example. In this case, the elemental structure is a one-bolometer Wheatstone bridge, and 
the Joule heating voltage is electronically subtracted using a dark current subtraction 


method. The thermal characteristics of the bolometer are given below. 





Performance Parameter Capability 


Pixel Dimensions 50 (um)* 
Thermal Conductance (G) <5x10° WK 
Thermal Capacity (H) <1x10° JK 
Thermal Time Constant (t) 20 msec 

Oo > 02 
Resistance (Ro) 50 kQ 


Table III.2 Raytheon IR imaging rifle-sight bolometer characteristics. 


Given such a bolometer, we would like to determine the advantage gained in using a 
compensating bolometer in the Wheatstone bridge configuration. From equation (II.15) 
we see that the saturation voltage from a one-bolometer Wheatstone bridge is given by 


Via P 
Ll cueaiee : 4 G . 


However, we must consider the saturation voltage for a Wheatstone bridge with an 





(111.8) 


additional compensating bolometer, whose thermal conductance is different from the 
sensor bolometer by definition. This calculation would be analytically impossible were it 
not for the principle of superposition (Boas, 1983), which allows us to analyze the 


Wheatstone bridge twice, once each with a different bolometer, and then add the results. 


The principle of superposition is valid for any linear network containing resistors and 
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voltage sources. Using this method in the situation where the compensating bolometer is 
identical to the sensor bolometer with the exception of thermal conductance, we retrace 


the steps which lead to equation (II.15) and arrive at 


. ahd) aL9) 


saturation A 
A convenient figure of merit for the detector would be its signal response at 4 nW of 
incident IR flux. For the Raytheon bolometer circuit the signal would be 1.2 mV at Vp=3 
V. For a reference bolometer with thermal conductance G; one order of magnitude 
greater (5x10°’ W/K), the signal would be 1.08 mV. There has been a 10% reduction in 
signal strength, and while the signal is certainly strong enough to manipulate, any 
reduction in signal strength is not desirable. However, the following graph illustrates the 
gain in integration time, assuming the signal will be integrated until the J oule heating 


matches the signal strength. 
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Figure III.13 | Comparison of integration times with and without a compensating 
bolometer. 
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Clearly the benefit in integration time far outweighs this loss of signal strength, when one 
considers the effect on noise that an integration time of approximately two microseconds 
versus approximately ninety microseconds has. Following equation (III.7), the one 
bolometer circuit has a bandwidth of 5x10° Hz while the compensated circuit has a 
bandwidth of 1.1x10* Hz. It follows that the noise in the compensated bolometer is one- 


and-a-half orders of magnitude less than in the uncompensated circuit. 


For imaging devices, this kind of noise reduction greatly enhances image quality. 
Video imagery is typically around 30 frames per second (fps), and even with the 
extended integration time we are nowhere near this limit. Even with the 10% loss of 
signal, implementation of a compensating bolometer in the IR imaging rifle-sight would 
improve it dramatically. 

3. Further Possible Benefit — Cancellation by Integration 

The benefit of having developed a working model of the compensated 
Wheatstone bridge becomes manifest when one simply begins to tinker with the design of 
the compensating bolometer. Instead of spending valuable resources in the design, 
manufacture and testing of physical bolometers, concepts can be modeled with ease. 
Such simulations may hold'the key to the actual elimination of the Joule heating output 
voltage. In adjusting the thermal capacities and conductivities somewhat haphazardly, it 
was found that a Joule heating curve could be generated which first went negative, and 
then curved back positive before it exploded into the normal exponential growth curve. 
This was intriguing because it appeared to offer a means of canceling out the Joule 
heating voltage; simply integrating for a period long enough to cancel the negative 


voltage area with an equivalent positive voltage area would eliminate the Joule heating 
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voltage. Unfortunately, it was noted that adjusting the thermal mass of the reference 
bolometer would be very impractical in the micro-machining process, and therefore 


prohibitively expensive. 


Fortunately, from the analytic description of the bolometer’s early time response 
another option for achieving the same behavior appears: 


3 
Va 


V, =—t———t. (11.19) 
16-R,-H 


Equation (11.19) indicates that output voltage is equally dependent on the balancing 
resistor values as the thermal mass. This leads to the idea that the same advantageous 
output voltage curve can be generated by manipulating the value of Ro on one side of the 


Wheatstone bridge. The figure below illustrates the idea. 


V 





R,(1+aAT) R,(1+aAT) 


F igure 111.14 Unbalanced Wheatstone bridge 


The design shown in figure III.14 was implemented in PSPICE, with the same 
Raytheon bolometer as described by table III.2. This time, instead of precisely balancing 
R; and R> to equal each other, the value of R2 was adjusted downward, resulting in the 


following set of output curves. 
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Figure IfI:15 PSPICE simulated output voltage versus time for integration by 
cancellation design. 


As hoped for, adjusting R2z lower than R; caused the curve to dip negative before 
launching exponentially in the positive direction. To integrate the Joule heating voltage, 
an integrator circuit, such as a Miller integrator, would have to be integrated into the pre- 
amplifier circuitry. It would therefore be important to keep the negative voltage dip 
roughly equivalent to the signal voltage, in order to keep the integrator’s capacitor small 
enough for an appreciable output signal. R2 meets this requirement, while providing for 
an integration time of approximately 450 microseconds. The Joule heating voltage could 
be entirely cancelled out, leaving behind only the signal voltage. Not only is the sensor 
greatly enhanced by removal of the Joule heating voltage, but the integration time has 


been extended, further improving the noise reduction by integration, and hence the 
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signal-to-noise ration (SNR). For example, the predominate source of noise, Johnson 


noise, may be calculated as follows. 


Using the principle of superposition, we can derive the noise generated by each of 
the four resistors individually, and then add the noise figures in quadrature, as should be 
done with uncorrelated noise sources. Assuming all noiseless resistors with one noise 


voltage source for the lower right resistor, the circuit appears below. 





Figure III.16 Wheatstone bridge with voltage noise source for lower right resistor 


The development is similar to that for equation (II.11), but this time V2 is different, with 





V,=-b—— (11.10) 


instead of V, = Ys oo where V, is the noise voltage due to one resistor. It then follows 


that 
AV cise =|V, V,| a rs mf 
2 2 
AV noise x “2 : (It. 1 1) 

















For the total noise voltage, the individual noise voltages from each resistor are added in 


quadrature. In this case, the four individual noise voltages are equivalent, so 


or 
Vwose =, - (11.12) 
In dealing with Johnson noise, 
Voie = V4kTRAS . (HII.13) 


We may now evaluate the benefit of the longer integration time provided by the 
cancellation by integration technique. From figure III.9 we conservatively estimate an 
integration time of t = 2 us, which also leaves behind 0.6 mV of Joule heating voltage. 
Using this integration time, R= 50 kQ and 7=294 K, equation (III.13) yields 


Vvorse = 14.25 wv. (111.14) 


Comparing this with cancellation by integration using R2 = 47kQ, we can estimate t = 
425 us from figure II.11. This time equation (II.13) yields 


Vvorse = 0.98 wV (111.15) 


In essence, these figures represent the minimum detectable signal threshold, because the 
sampling time is still short enough that A, noise is negligible compared to Johnson 
noise. In addition, the normally remaining Joule heating voltage can be completely 


eliminated, further reducing the inaccuracies that accompany electronic subtraction of the 


Joule heating voltage. 


This appears to be a very promising technique, because, while altering an 
individual bolometer’s thermal mass (H) is very difficult (requiring a three-dimensional 


architectural change), altering its thermal conductance (G) and baseline resistance is 
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simply a matter of changing the two-dimensional architecture of the reference bolometer 


in the micro-machining process. 
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IV. CONCLUSION 
By developing a working model on PSPICE, we have demonstrated the ability to 


explore more deeply the complex interaction between the thermal characteristics of the 
paired bolometers in the compensating Wheatstone bridge. Understanding the complex 
thermal interactions is only possible because the computer model numerically accounts 
for the higher-order terms which are neglected by necessity in analytic treatments. In 
practical terms, the model allows the researcher to tinker and to explore these interactions 


without incurring enormous costs or wasting valuable research time. 


One example of the benefit of this kind of free tinkering was arrived at very 
quickly, with the discovery of cancellation by integration while using an unbalanced 
bridge. Previous thought had been entirely directed at the balanced bridge. Through the 
cancellation by integration technique developed in section III.C.3, we were able to realize 
a dramatic increase in integration time, which equated to a minimum detectable signal 


more than an order of magnitude smaller than without the technique. 


A great deal of information awaits those who wish to more thoroughly explore the 
Wheatstone bridge arrangement using this model. With five independent variables all 
affecting the output voltage, there remain many unusual configurations worthy of study. 
One obviously worthwhile pursuit would be to study the sensitivity to manufacturing 
tolerances, in order to determine how tightly a particular design must be manufactured to 
deliver a tolerable performance. Another beneficial pursuit would be that of 
understanding the physical limits of a microbolometer based on the architectural design 
limits of the micromachining process (i.e., how narrow a bridge can be made and what 


this means in terms of thermal conductance limits). 


a1 
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APPENDIX A: THE BOLOMETER ARRAY 


A bolometer array can be designed with the reference bolometer outside of the 
pixel array. Using vertical and horizontal shift registers, each individual bolometer can 
be read sequentially and incorporated into the Wheatstone bridge structure in the lower 


right corner. 
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Ee att > 
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Figure A.1 |= Microbolometer Focal Plane Array (FPA) architecture. 


Vertical Shift Register 
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